Superfluid transition in quasi2D Fermi gases 
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We show that atomic Fermi gases in quasi2D geometries are promising for achieving superfluidity. In 
the regime of BCS pairing for weak attraction, we calculate the critical temperature T c and analyze 
possibilities of increasing the ratio of T c to the Fermi energy. In the opposite limit, where a strong 
coupling leads to the formation of weakly bound quasi2D dimers, we find that their Bose-Einstein 
condensate will be stable on a long time scale. 
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Recent progress in trapping and cooling of Fermi iso- 
topes of K (Hi and Li has shown the ability 
to go far below the temperature of quantum degeneracy 
and to manipulate independently the trapping geometry, 
density, temperature and interparticle interaction. The 
Duke experiment M presents intriguing results on the 
possibility of achieving a superfluid phase transition in 
the two-component Fermi gas of 6 Li. 

Two-dimensional Fermi gases have striking features 
not encountered in 3D. In the superfluid state, thermal 
fluctuations of the phase of the order parameter strongly 
modify the phase coherence properties. The interaction 
strength depends logarithmically on the relative energy 
of the colliding atoms. For degenerate Fermi gases this 
energy is of the order of the Fermi energy £p which is 
proportional to the 2D density n. Accordingly, the expo- 
nential dependence of the BCS transition temperature on 
the interaction strength transforms into a power law de- 
pendence on the density: T c oc v}l 2 [OUT!]. This suggests 



a unique possibility to cross the critical point by adiabat- 
ically expanding a degenerate Fermi gas. Since the ratio 
T/ep remains unchanged, the temperature scales as n 
and decreases with density faster than T c . 

Experimentally it is possible to achieve the quasi2D 
regime by confining the atoms in one direction so tightly 
that the corresponding level spacing exceeds the Fermi 
energy. Under this condition the degenerate Fermi gas is 
kinematically two-dimensional. Thus far, this regime has 
been reached for Cs atoms pj9|, |lo| and for Bose-Einstein 
condensates of Na and Rb |12| . 

In the quasi2D regime the mean-field interaction be- 
tween particles exhibits a similar logarithmic dependence 
on the particle energy as in the purely 2D case |l5|] . The 
amplitude of the s-wave scattering turns out to be sen- 
sitive to the strength of the tight confinement jl5| . This 
opens new handles on manipulations of the interparticle 
interaction and superfluid pairing. 

In this Rapid Communication we show that atomic 
Fermi gases in quasi2D geometries can become strong 



competitors of 3D gases in achieving superfluidity. The 
ability to increase the interparticle interaction by tun- 
ing the trap frequencies gives an opportunity to realize 
a transition from the standard BCS pairing in the case 
of weak attraction to the limit of strong interactions and 
pairing in coordinate space. In the latter case one eventu- 
ally gets a dilute system of weakly bound quasi2D dimers 
of fermionic atoms, which can undergo Bose-Einstein 
condensation. For the BCS case, we calculate the critical 
temperature T c to second order in perturbation theory 
and discuss possibilities of increasing the ratio T c /ef- In 
the other extreme, we find that the interaction between 
the quasi2D dimers is repulsive, and their collisional re- 
laxation and decay are strongly suppressed. This allows 
us to conclude that BEC of these composite bosons will 
be stable on a long time scale. 

We consider an ultracold two-component Fermi gas in 
the quasi2D regime and confine ourselves to the s-wave 
interaction and superfluid pairing between atoms of dif- 
ferent components. We assume that the characteristic 
radius of the interaction potential is much smaller than 
the harmonic oscillator length in the tightly confined di- 
rection, la = (h/mtoo) 1 ' 2 , where m is the atom mass, 
and wo is the confinement frequency. Then the interac- 
tion problem involves two length scales: lo and the 3D 
scattering length a. For a < and \a\ <SC lo, there is a 
peculiar quasi2D weakly bound s-state of two particles, 
with the binding energy 15 



£o = 0.915(7kjo/tt) ex P (— V 27rZ /|a|) <C Tiujq 



(1) 



In this case the coupling constant for the intercompo- 
nent interaction takes the form g = (4-irh 2 /m) ln _1 (eo/e), 
where the relative collision energy e is assumed to be ei- 
ther much smaller or much larger than £o (see |15 and 
refs. therein). As in degenerate Fermi gases one has 
e ~ Ef, the interaction is attractive (g < 0) if the den- 
sity is sufficiently high and one satisfies the inequality 



E /e F < 1. 



(2) 
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Thus, the inequality is the necessary condition for 
the BCS pairing. For finding the critical temperature T c 
below which the formation of Cooper pairs becomes fa- 
vorable, we go beyond the simple BCS approach and pro- 
ceed along the lines of the theory developed by Gor'kov 
and Melik-Barkhudarov for the 3D case |l6| . 

The critical temperature T c is determined as the high- 
est temperature for which the linearized equation for the 
order parameter (gap) A = (g^^) has a nontrivial solu- 
tion |l7|]. Assuming that the quasi2D gas is uniform in 
two in-plane directions, the gap equation in the momen- 
tum space takes the 2D form 



A(q) q,*) 



K(q') 



1 



+ 5V( q ,q')K( q ')\A( q >)- d2 ' 1 



(27T) 



(3) 



where K{q) 



(l/2£(g))tanh(£(g)/2T), £(<?) 



H q /2m — fx, and fi « ef = ffTi n/m is the chemical po- 
tential. The first term in the rhs of Eq.(j3|) results from 
the direct interaction between particles, and we renor- 
malized the interaction potential in terms of the scatter- 
ing amplitude / (vertex function). The latter is a solu- 
tion of the quasi2D scattering problem. The parameter 
z has a meaning of the total energy of colliding particles 
in their center of mass reference frame. It is of the or- 
der of Ep and drops out of the final answer. The term 
5V(q, q') describes the modification of the interparticle 
interaction due to the presence of other particles (many- 
body effects). The leading contributions to this term 
are second order in the scattering amplitude. They are 
shown in Fig^ and correspond to an indirect interaction 
between two particles when one of them interacts with a 
particle-hole pair virtually created from the ground state 
(filled Fermi sea) by the other particle. These second or- 
der contributions are important for the absolute value of 
the critical temperature (preexponential factor in the 3D 
case), whereas higher order terms involving more inter- 
action events can be neglected. 
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FIG. 1. The leading contributions to 5l / (q, q'). 



The amplitude / is independent of the momenta q, q' 
and we will use the 2D relation (see Ea]) 



f = A-K\zT 1 {E /{-z)}. 



(4) 



The last term in the rhs of Eq. (||) is a small correction, 
since SV ~ / 2 . Accordingly, the quantity |/| represents 
a small parameter of the theory. In the quasi2D regime 
(z <C hu>o) the motion of particles in the tightly con- 
fined direction provides a correction to Eq.(^), which is 
~ (z/fnuo)/ 2 Jl5fl . It is much smaller than SV and will 
be omitted. 

Equations (||) and (Q) show that the momentum de- 
pendence of the order parameter appears only due to 
the second order term that contains many-body contri- 
butions to the interparticle interaction SV. The latter 
is a function of p — |q + q'| and rapidly decays for 
p > 2qp, where qp — y/2~mJl/Ti is the Fermi momen- 
tum. For p < 2qp the quantity SV is almost constant. 
Therefore, one has A(q') ~ A(qp) in a wide momentum 
range near the Fermi surface. Then, for q = qp a direct 
integration of Eq. (0) yields 



A(g F ) = -^A( 9F )ln 

TYl 

- SV{qF,q F )—^ 
znn 



exp(2 7 ) 



7T 2 T 2 

ln(c|) A(q F ), 



where 7 » 0.5772 is the Euler constant, and C is a nu- 
merical factor determined by the momentum dependence 
of A and SV. The calculation of 5V(qF,qp) is straight- 
forward and gives SV(q F ,q F ) = {h 2 /2Trm)f 2 (2fj,). Then, 
using Eq.(|J) we obtain the critical temperature 

T c « (2 A1 / 7 r)exp (7 - 1 - |2 7 rRc/- 1 (2 A1 )|) . 

The exponent in this equation should be large and the 
quantity Re/ _1 (2/i) should be negative. As the chemi- 
cal potential is fj, w Ef, from Eq.(Q) one sees that these 
requirements are reached under the condition (|^). Using 
Eqs. (Q) and (ji|) the critical temperature takes the form 



7v / 2e u 7F 



0.16VeF^oexp (r\J^T^j\ < £ f- (5) 



The relative correction to this result is of the order of 

l/|ln(e /£ F )| « 1. 

Note that Eq. (|B|) predicts by a factor of e smaller 
value for the critical temperature than a simple BCS cal- 
culation |n| . This means that the attractive interaction 
between particles becomes weaker once we take into ac- 
count the polarization of the medium. 

The ratio T c /s p is not necessarily very small. For ex- 
ample, using Feshbach resonances the scattering length 
is tunable over a wide interval of negative values ■ 
Keeping the exponential term equal to 0.05 in Eq.(pp, 
with u>q ~ 100 kHz we obtain T c /sf ~ 0.1 for 2D densi- 
ties n ~ 10 9 cm" 2 (T c ~ 40 nK). 
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As in the purely 2D case JL3|J14| ], the transition tem- 
perature T c oc n 1 / 2 and the ratio T c /£p increases with 
decreasing density as rT 1 / 2 . This is a striking differ- 
ence from the 3D case, where this ratio decreases expo- 
nentially with density. In the presence of the in-planc 
confinement, one can approach the BCS transition in 
a degenerate Fermi gas by adiabatically expanding the 
quasi2D trap in the in-plane direction(s). As the degen- 
eracy parameter T/sp is conserved in the course of the 
adiabatic expansion, the ratio T/T c will decrease as n 1 ^ 2 . 
Equations ([!]) and (||) also show that one can increase Eq 
and T c /e by tuning \a\ to larger values or by making the 
tight confinement stronger and thus decreasing Iq. 

What happens if eo and Ef become comparable with 
each other, i.e. one reaches the quasi2D resonance for 
two-body collisions? Then Eq.(|J) leads to T c ~ £j? and 
is no longer valid. In fact, for £q > ef the forma- 
tion of bound quasi2D dimers of distinguishable fermions 
becomes energetically favorable and one encounters the 
problem of Bose-Einstein condensation of these bosonic 
molecules. Thus, an increase of the ratio Eq/ef from 
small to large values is expected to provide a transfor- 
mation of the BCS pairing to molecular BEC. This type 
of crossover has been discussed in literature in the con- 
text of superconductivity pO| , pT| , p2p^,p^ ] and in relation 
to superfluidity in 2D films of 3 He [ |l3|j24[ ]. The idea of us- 
ing a Feshbach resonance for achieving a superfluid tran- 
sition in the BCS-BEC crossover regime in ultracold 3D 
Fermi gases has been proposed in refs. p^j2^ |. 

We will not consider the crossover regime and con- 
fine ourselves to the limiting case of molecular BEC 
(eo ^> £f)- A subtle question is related to the stabil- 
ity of the expected condensate, which depends on the 
interaction between the molecules. For the repulsive in- 
teraction one will have a stable molecular BEC, and the 
attractive interaction should cause a collapse. 

The molecule-molecule scattering is a 4-body problem 
described by the Schrodinger equation 



- v 

m 



x-V 2 n + U( ri ) + U(r 2 ) 
2m 



E 



*(r 1; r 2 ,R) = 0. 



(6) 



Here i*i is the distance between two given distinguishable 
fermions, r 2 is the distance between the other two, R is 
the distance between the centers of mass of these pairs, 
and U is the interatomic potential. The total energy is 
E = —2eq+£, with e being the relative molecule-molecule 
kinetic energy. 

The interaction between molecules is present only 
at intermolecular distances of the order of or smaller 
than the size of a molecule = %/ ' ^/ttieq. Therefore, 
at energies e <C £o the scattering between molecules 
is dominated by the s-wave channel and can be ana- 



lyzed on the basis of the solution of Eq.(||) for e = 
0. For large R the corresponding wavefunction is 
*(ri,r 2 ,R) w Ko(ri/d*)K {r 2 /d*)bx(aR/di,), where 
the decaying Bessel function Ko(ri t 2/d*) represents the 
2-body bound state. The parameter a is a universal 
constant which can be found by matching the quan- 
tity ln(ai?/d*) with the solution of Eq.(^) at short dis- 
tances. Finally, matching ln(ai?/d*) with the wavefunc- 
tion of free relative motion of two molecules at distances 
d* <C R <C A £ , where A e = %/ %Jvtie is their de Broglie 
wavelength, we obtain the coupling constant (scattering 
amplitude) for the interaction between molecules: 

g m = {2nh 2 /m) In" 1 (2a 2 e" 27 £ /e) > 0; £ < e . (7) 

A precise value of a is not important as it gives rise 
to higher order corrections in Eq.(Q). However, in or- 
der to make sure that this constant is neither anoma- 
lously large nor anomalously small we have integrated 
Eq.(^|) numerically. For this purpose, it is convenient to 
transform Eq. (g) into an integral equation for a function 
which depends only on three independent coordinates. 
This has been done by using the method of ref. p7| . 
Our calculations lead to a 1.6. They show the ab- 
sence of 4-body weakly bound states and confirm an intu- 
itive picture that the interaction between two molecules 
can be qualitatively represented by means of a purely 
repulsive potential with the range ~ . For the interac- 
tion between Bose-condensed dimers, in Eq.(^) one has 
£ = 2n m g m <C £o, where n m is the density of the dimers 
(see |D| and refs. therein). We thus conclude that a 
Bose condensate of these weakly bound dimers is stable 
with respect to collapse. 

The 2D gas of bosons becomes Bose-condensed below 
the Kosterlitz-Thouless transition temperature Tkt pS|] 
which depends on the interaction between particles. Ac- 
cording to the recent quantum Monte Carlo simulations 
p9| , for the 2D gas with the coupling constant (Q) the 
Kosterlitz-Thouless temperature is given by 

T KT = (Trh 2 n m /m) In" 1 [( V /4tt) In (l/n m d 2 )] , (8) 

where the numerical factor r\ w 380. For ep <C £o, 
the density of dimers n m rs n/2 and the parameter 
{l/n m d 2 ) ps 2-keo/ef- Then, {orso/sp = 10, Eq.(||) gives 
Tkt/ef 0.1 and at densities 10 8 cm -2 the transition 
temperature in the case of 6 Li is Tkt ~ 30 nK. 

The weakly bound dimers that we are considering are 
molecules in the highest rovibrational state and they can 
undergo collisional relaxation and decay. The relaxation 
process occurs in pair dimer-dimer or dimer-atom col- 
lisions. It produces diatomic molecules in deep bound 
states and is accompanied by a release of the kinetic en- 
ergy. The size of these deeply bound molecules is of the 
order of the characteristic radius of the interatomic po- 
tential R e Iq, and their internal properties are not 
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influenced by the tight confinement. Therefore, the re- 
laxation can be treated as a 3D process and it requires 
the presence of at least three fermionic atoms at dis- 
tances ~ R e between them. Since at least two of them 
are identical, the relaxation probability acquires a small 
factor (kR e ) 2 compared to the case of bosons, where 
k ~ 1/d* = e$ / Tiujq / Iq is a characteristic momentum of 
atoms. The 3D density of atoms in the quasi2D geometry 
is ~ ti/Iq. Thus, qualitatively, the inverse relaxation time 
can be written as r~j ~ a le in(R e /lo) 2 (£o/^o)/lo, where 
a rc i is the relaxation rate constant for the highest rovi- 
brational states of 3D molecules of two bosonic atoms. 
We estimate r ro i keeping in mind the recent measure- 
ments for Rb2 molecules ||(]] which give a Te \ f=s 3 x 10 -11 
cm 3 /s. For l in the interval from 10 -5 to 10 -4 cm, the 
suppression factor (R e /l ) 2 (so/huJo) ranges from 10~ 3 to 
10 -5 and at 2D densities n ~ 10 s cm' 2 we find the re- 
laxation time T re i of the order of a second or larger. 

Dimcr-dimer pair collisions can lead to the formation 
of bound trimers, accompanied by a release of one of the 
atoms. The formation of deeply bound trimer states will 
be suppressed at least in the same way as the relaxation 
process discussed above. Therefore, it is important that 
there are no weakly bound trimers in (quasi) 2D. Just as 
m 3D this can be established by using the zero-range 
approach [R e — > 0). We have performed this analysis 
along the lines of the 3D work p7j . Qualitatively, the 
symmetry of the 3-fermion system containing two identi- 
cal fermions provides a strong centrifugal repulsion that 
does not allow the presence of 3-body bound states. This 
is in contrast to 2D bosons where one has two fully sym- 
metric trimer bound states |32] ]. 

Thus, the life-time of quasi2D dimers of fermionic 
atoms is rather long and one easily estimates that it 
greatly exceeds the characteristic time of elastic colli- 
sions. One can even think of achieving BEC in the ini- 
tially non-condensed gas of dimers produced out of a non- 
superfluid atomic Fermi gas under a decrease of n or Iq. 

In conclusion, we have found the temperature of super- 
fluid phase transition in two-component quasi2D Fermi 
gases. Our results are promising for achieving this tran- 
sition in both the regime of BCS pairing and the regime 
of BEC of weakly bound dimers. 
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